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Abstract. For triangulated surfaces, we introduce the combinatorial Calabi flow which 
is an analogue of smooth Calabi flow. We prove that the solution of combinatorial Calabi 
flow exists for all time. Moreover, the solution converges if and only if Thurston's circle 
packing exists. As a consequence, combinatorial Calabi flow provides a new algorithm 
to find circle packings with prescribed curvatures. The proofs rely on careful analysis of 
combinatorial Calabi energy, combinatorial Ricci potential and discrete dual-Laplacians. 



1. Introduction 

In seeking constant curvature metric, E. Calabi studied the variational problem of mini- 
mizing the so-called "Calabi energy" in any fixed cohomology class of Kahler metrics and 
proposed the Calabi flow ([2], [3]). Hamilton introduced the Ricci flow technique (|19j). 
which has been used to solve the Poincare conjecture. As to dimension two, i.e., smooth 
surface case, people proved that both Calabi flow and the normalized Ricci flow exist for 
all time t > 0. Moreover, these two flows converge to constant scalar curvature metric 
(see @], 0, @, 0, 1, [21!, and [25]). 

Given a triangulated surface, Thurston introduced the circle packing metric, which is a 
type of piecewise flat cone metric with singularities at the vertices. Thurston found that 
there are combinatorial obstructions for the existence of constant combinatorial curvature 
circle packing metric (see section 13.7 in |26j). In [8j, Bennett Chow and Feng Luo defined 
combinatorial Ricci flow, which is the analogue of Hamilton's Ricci flow. Bennett Chow 
and Feng Luo proved that the combinatorial Ricci flow exists for all time and converges 
exponentially fast to Thurston's circle packing on surfaces. They reproved the equivalence 



between Thurston's combinatorial condition (see (1.3) in |8j, or (1.3) in this paper) and 
the existence of constant curvature circle packing metric. 

Inspired by the work of Bennett Chow and Feng Luo in [8], we consider in this paper 
the 2-dimensional combinatorial Calabi flow equation, which is the negative gradient flow 
of combinatorial Calabi energy. We interpret the Jacobian of the curvature map as a 
type of discrete Laplace operator, which comes from the dual structure of circle patterns. 
We evaluate an uniform bound for discrete dual-Laplacians, which implies the long time 
existence of the solutions of combinatorial Calabi flow. We prove that the combinatorial 
Calabi flow converges exponentially fast if and only if Thurston's combinatorial conditions 
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Figure 1. circle packing metric 



Figure 2. two adjacent triangles 



are satisfied. The combinatorial Calabi flow tends to find Thurston's circle patterns auto- 
matically. As a consequence, we can design algorithm to compute circle packing metrics 
with prescribed combinatorial curvature. In fact, any algorithm minimizing combinatorial 
Calabi energy or combinatorial Ricci potential can achieve this goal. 

1.1. Circle packing metric. Let's consider a closed surface X. A triangulation T = 
(V,E,F) has a collection of vertices (denoted V), edges (denoted E), and faces (denoted 
F). A positive function r : V — > (0, +oo) defined on the vertices is called a circle packing 
metric and a function $ : E — )■ [0, 7r/2] is called a weight on the triangulation. Throughout 
this paper, a function defined on vertices is regarded as a iV-dimensional column vector 
and N = V^ 1 is used to denote the number of vertices. Moreover, all vertices, marked by 
ui, ujv, are supposed to be ordered one by one and we often write i instead of Vi. Thus 
we may think of circle packing metrics as points in R> , N times of Cartesian product of 
(0, oo). A triangulated surface with wight $ is often denoted as (X, T, <I>). 

For fixed (X, T, $), every circle packing metric r determines a piecewise linear metric 
on X by attaching each edge to a length 



This length structure makes each face in F isometric to an Euclidean triangle, more 
specifically, we can realize each face {i, j, k} E F as an Euclidean triangle with edge 
lengths kj,ljk,hi because the three positive numbers kj,ljk,hi (derived from r%, rj, r&) 
must satisfy triangle inequalities ([2S]) Lemma 13.7.2). Furthermore, the triangulated 
surface (X, T) is composed by gluing many Euclidean triangles coherently. 

1.2. Combinatorial curvature and constant curvature metric. Given (X, T, <3?) and 
a circle packing metric r, all inner angles of the triangles are determined by r%, ■ ■ ■ ,r^. 
Denote 9{ k as the inner angle at vertex i in the triangle {i, j, k} G F, then the well known 
combinatorial (or "discrete") Gauss curvature Ki at vertex i is 





(1.1) 



{i,j,k}eF 
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where the sum is taken over each triangle having i as one of its vertices. Notice that 9j k 
can be calculated by cosine law, thus 9j k and Ki are elementary functions of the circle 
packing metric r. 

For every circle packing metric r on (X, T, $), we have the combinatorial Gauss-Bonnet 
formula ([8]) 

JV 

Y,Ki = 2n X (X). (1.2) 
i=i 



Notice that, above combinatorial Gauss-Bonnet formula (1.2) is still valid for all piecewise 
linear surface (X, T, I), where I : E — > (0, +oo) is a piecewise linear metric, that is, the 
length structure I makes each face in F isometric to an Euclidean triangle. The average 
combinatorial curvature is k av = 2ttx(X)/N . It will not change as the metric (circle 
packing metric or piecewise linear metric) varies, and it's a invariant that depends only 
on the topological (the Euler characteristic number x(X)) and combinatorial (N = V^ 1 ) 
information of (X, T) . 

Finding good metrics is always a central topic in geometry and topology. The constant 
curvature circle packing metric (denoted as r av ), a metric that determines constant combi- 
natorial curvature K av = K(r av ) = k av (l, • • • , 1) T , seems a good candidate for privileged 
metric. Thurston first studied this class of metrics, and found that there are combinatorial 
obstructions for the existence of constant curvature metric (|26j). We quote Thurston's 
combinatorial conditions here. 

Lemma 1.1. Given (X, T, where X, T and $ are defined as before. Then there 
exists constant combinatorial curvature circle packing metric if and only if the following 
combinatorial conditions are satisfied 

2tt X (X)^- >- Yl (T-*(e))+27r X (F / ), VI : <t> + 1 C V, (1.3) 

where Fj is the subcomplex whose vertices are in /, Lk(I) is the set of pairs (e, v) of an 
edge e and a vertex v so that (1) the end points of e are not in I and (2) the vertex v is 
in / and (3) e and v form a face in T. 

In [8] , Bennett Chow and Feng Luo introduced the combinatorial Ricci flow which is an 
analogue of Hamilton's Ricci flow on surfaces ([?]). They proved that the solutions of the 
normalized combinatorial Ricci flow converges exponentially fast to Thurston's constant 
curvature metric r av . Combinatorial Ricci flow method is so powerful that, on one hand, 
the flow provides a new proof of Thurston's existence of constant circle packing metric 
theorem, on the other hand, the flow suggests a new algorithm to find circle packing 
metrics with prescribed curvatures. 

In this paper, we introduce the combinatorial Calabi flow, which is the negative gradient 
flow of combinatorial Calabi energy. Combinatorial Calabi flow is an analogue of smooth 
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Calabi flow on surface. We prove that the solution of combinatorial Calabi flow exists for 
t G [0, +00) by a careful estimation of discrete dual-Laplacian. Moreover, the solution of 
combinatorial Calabi flow converges to constant curvature metric r av if and only if r av 
exists. 



2. The 2-Dimensional Combinatorial Calabi Flow 

2.1. Definition of combinatorial Calabi flow. For smooth surface, the well known 
Calabi flow is 3| = AKg, where K is Gauss curvature. Obviously, smooth Laplace- 
Beltrami operator is very important for smooth Calabi flow equation. Before giving the 
definition of combinatorial Calabi flow, we first need to define combinatorial Laplace op- 
erator, which is an analogue of smooth Laplace-Beltrami operator. 

Set U{ = lnrj, i = 1, • • • ,N. Then the coordinate transformation u = u{r) maps 
r G R^g to u G H homeomorphically (denote the inverse coordinate transformation as 
Exp : R N ™ N 

Laplacian to be the Jacobian of the curvature map K 



H> , u 1 y t for latter use). Similar to [11] and [12] , we interpret the discrete 



K(u). 

Definition 2.1. Given (X, T, <3?), where X, T and $ is defined as before. The discrete 
dual-Laplace operator "A" is defined to be A = — L T , where 



L = (L. 



ij)NxN 



d{K u ...,K N ) 
d(m, ...,u N ) 



/ dKi 
' dui 



OKn 



8K1 

dupf 



dK N 
du N 



\ 



Both A and L acts on functions / (defined on vertices, hence is a column vector) by 
matrix multiplication, i.e. 

N 

A/i = (A/)j = —(L T f)i 



Discrete Laplacian defined here is exactly the famous discrete dual-Laplacian ([13J [15] ). 
which will be explained carefully in Appendix |A1| 

Definition 2.2. Given data (X, T, $), where X, T and <3? is defined as before. The 
combinatorial Calabi flow is 



dt 



TiAKi 



with initial circle packing metric r(0) G I&>0] or equivalently, 

~dt 

with n(0) G R N . 



(2.1) 



(2.2) 
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It's more neatly to write combinatorial Calabi flow (2.1) and (2.2) in matrix forms. 
Denote R = diag{ri, 

(2.3) 



, r^v}, then flow (2.1) turns to 
r 



RAK 



RL T K, 



where the symbol " • " represents the time derivative. The flow (2.2) turns to 

u = AK = —L T K. 



(2.4) 



Notice that, equation (2.4) is an autonomous ODE system. 



Remark 1. Using the matrix language, Bennett Chow and Feng Luo's combinatorial 
Ricci flow ([8]) is u = — K. The normalized combinatorial Ricci flow is u = K av — K. 

2.2. Main properties of combinatorial Calabi flow. 

Definition 2.3. For data (X, T, <£) with a circle packing metric r, the combinatorial 
Calabi energy is 



N 



C(r) 



\K — K n 



(2.5) 



i=l 



Consider the combinatorial Calabi energy C as a function of u, we obtain 



/ 9C_ 
1 d Ul 



8C 



/ 9Ki 
1 d Ul 



J 



dKi 



Proposition 2.4. Combinatorial Calabi flow (2 



8K 



\ 



du N / 



2L T K. 



2) (or (2.4)) is the negative gradient flow 



of combinatorial Calabi energy. Moreover, the calabi energy (2.5) is descending along this 
flow. 

Proof: u(t) = —L T K = -\V U C. Moreover, C'(t) = -i||V u C|| 2 < 0. 

Q.E.D. 

Because Ki and AKi are elementary functions of ri, rjv, the local existence of com- 
binatorial Calabi flow (2.1) follows from Picard's existence and uniqueness theorem in the 
standard ODE theory. By a careful estimation of discrete dual-Laplacian and combinato- 
rial Ricci potential, the long time existence property is proved in section[3j the convergence 
property is proved in section [4} The main results in this paper is stated as follows 

Theorem 2.5. Fix (X, T, <3?), where X is a closed surface, T is a triangulation, $ S [0, n/2] 
is a weight. For any initial circle packing metric r(0) G B,>q> the solution of combinatorial 
Calabi flow (2.1) exists for t G [0, +oo). Moreover, the following four statements are 
mutually equivalent 



(1): The solution of combinatorial Calabi flow (2.1) converges. 
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(2) 
(3) 
(4) 



The solution of combinatorial Ricci flow u = K av — K converges. 
There exists constant curvature circle packing metric r a 



av ■ 



Furthermore, if any of above four statements is true, the solution of combinatorial Calabi 
flow (2.1) converges exponentially fast to constant curvature circle packing metric v av . 

3. Long Time Existence 

3.1. Discrete Laplacian for circle packing metrics. Any circle packing metric r de- 
termines an intrinsic metric structure on fixed (X, T, <&) by Euclidean cosine law. The 
lengths lij, angles 9j and curvatures Ki are elementary functions of r = (n, r^) T . We 
write j ~ i in the following if the vertices i and j are adjacent. For any vertex i and any 
edge j ~ i, set 

Q(Qi k + o jl ) 

By = 1 ^ ' J r„ (3.1) 

then = Bjj, since -g^-fj = -g^-ri (see Lemma 2.3 in jS]). 
Proposition 3.1. For any 1 < i, j < N and i ~ j, we have 

< < 2^3. (3.2) 

89 jk i- 

Proof: We just need to prove < -gyrrj < v3. We defer the details to the appendix. 

Q.E.D. 

Proposition 3.2. For any 1 < i,j < N, 

B ik , j = i 



-B^ , j ~ i ( 3 - 3 ) 
, else. 



Proof: This can be proved by direct calculations and we omit the details. 



Q.E.D. 



Proposition 3.3. L is semi-positive definite, having rank N — 1. Moreover, the kernel of 
L is the span of the vector (1, • • • , 1) T . 

Proof: This is a direct consequence from Lemma 3.10 in [8j. Q.E.D. 

The differential form oj = Yld=ii^i ~ k av )dui is closed, for Lij = Lji. Thus the integral 
Iuq Yli=i {^i ~ k a v)dui makes good sense ([H]), where uq is an arbitrary point in M, N . This 
integral is crucial for the proof of our main theorem. For convenience, we call this integral 
combinatorial Ricci potential. 

For any smooth closed manifold (M,g) with Riemannian metric g, we know that 
J M Af = 0, where / is an arbitrary smooth function defined on M. For combinatorial 
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surface (X, T, <I>) with circle packing metric r, we have A-fQ = — (1, • • • , 1) T LK = 0. 

i=l 

Thus we obtain 

Proposition 3.4. As long as the combinatorial Calabi flow exists, both Y\iLi r i(t) = 

Ilili r i(°) anci Y^iLi n i(0 = 12i=i are constants. 

3.2. Long time existence of combinatorial Calabi flow. Notice that 

AK i = Y,B ij {K j -K i ). (3.4) 

Using the estimation of in (3.2), we obtain 

Theorem 3.5. Given X, T, <3?, where X, T and 3> are defined as before. For any initial 
circle packing metric r(0), the solution of the combinatorial Calabi flow (2.1) exists for all 
time t G [0, +oo). 

Proof: Let d% denote the degree (or say valence) at vertex Vi, that is the number of 
edges adjacent to u$. Suppose d = max(d\ , ... , djv), evidently we see (2 — d)ir < Ki < 2tt. 
From the evaluation of Bij in Lemma 1 we know, all |Ai£j| is uniformly bounded by a 
positive constant c = 2\/3 ■ N ■ d ■ tt which depend only on the triangulation. Then we 
obtain 

c e~ ct < n{t) < c e ct 

where cq = c(r(0)) is a positive constant depends only on the initial state of r{t). This 
implies that the combinatorial Calabi flow has a solution for all time t S [0, oo) for any 
initial value r(0) € E> . Q.E.D. 

Remark 2. The long time existence of r{t) can be deduced without using the uniform 
estimation of Bij in (3.2). See Remark [3j 

4. Convergence to Constant Curvature Metric 

4.1. Evolution of combinatorial curvature. We consider a general flow in the space 
of all circle packing metrics with form 

dti 

where ifi are smooth functions of n, ...,rjv. For this flow, we have 
Lemma 4.1. Under general flow, the discrete Gauss curvature evolves according to 

= -A(fi = (Ltp)i 

or K = Lip in matrix form. 
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Proof: We can calculate —h- and -gf- in components form to get the proof. The 
detail is tedious and is omitted. However, above Lemma is a direct conclusion of chain 
rule in matrix form 

■ d(K lr .- ,K N ) . 
K = — -u = Lip = -Aip. 

Q.E.D. 



As to combinatorial Calabi flow, select ipi = AKi (or tp = —LK), and use lemma 4.1 
we obtain 

Proposition 4.2. Under combinatorial Calabi flow, the discrete Gauss curvature evolves 
according to 



4.2. Global rigidity and admissible space of curvatures. For fixed data (X, T, $), 

>N 
L >0 



let IT : R^g — > M, N be the curvature map which maps r to K. For any positive constant c, 



denote 

N 

^c={r=(r 1 ,--. ,r N ) T e^ \l[r t = c}. ^ 

i=l 

Thurston proved that the map II ^> c is injective, i.e., the metric is determined by its 
curvature up to a scalar multiplication ([26], [5]). This phenomenon is called "global 
rigidity". Denote 

jeb B = {K = (K lr -- ,K N ) T eR N | ^^ = 2vr X (X)}. 

i=l 

For any nonempty proper subset / of vertices V, set 

N 

m = {(K = K 1 ,--- ,K N ) T €R N \ Y, K i>~ E (vr-$(e)) + 2vr X (F / )} 

Let 

f ^jf GJj n(n w f;), (4.3) 

Bennett Chow and Feng Luo proved that ^ = H (H>o) ([5])' so ' ^ ^ s exactly the space 
of all possible combinatorial curvatures. We call W "admissible space of curvatures" . 
Notice that <3f is a nonempty bounded open convex polyhedron in (iV — l)-dimensional 
hyperplane J^gb- Moreover, & is homeomorphic to R^ -1 . By the invariance of domain 
theorem in differential topology, the curvature map II realized a homeomorphism from & c 
to <3f. Denote 

TV 

W a = {u = (u 1 ,...,u N f €B. N \ £)«i = a}, (4.4) 

i=i 
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then <^ c = Exp(^ a ), where a = Inc. Denote II = II o Exp : R 



N 



R 



JV 



u 



(->■ K. Then 



we have homeomorphisms % a P ^ \ 

positive constant c, & = H (R>o) = ^ 



ri-^K=n(r) 



n(R 



JV\ 



> ^ and, moreover, for arbitrary 
- II where a = Inc. 



4.3. Convergence to constant curvature metric. The constant curvature K a 
k av (l, ••• ,1) T may not admissible, that is, K av may not belongs to <3f . K a 
equivalent to the existence of constant curvature metric r av , actually, Ti~ l (K av ) 
a class of constant curvature metrics which differ by scalar multiplications. 



G 9 is 
is 



cR^o 



Along the combinatorial Calabi flow (2.1), from Proposition 3.4 we know that {r(£)} C 
where c = ]X=i ?™i(0). Now we are at the stage of proving convergence results. We 
state our main results as follows: 

Theorem 4.3. Fix (X, T, $), where X is a closed surface, T is a triangulation, $ S [0, n/2] 
is a weight. For any initial circle packing metric r(0) £ R>o> the solution of combinatorial 



Calabi flow (2.1) exists for t E [0, +oo). Moreover, the following four statements are 



mutually equivalent 



(1) 
(2) 
(3) 
(4) 



The solution of combinatorial Calabi flow (2.1) converges. 
The solution of combinatorial Ricci flow u 



K av — K converges. 
There exists constant curvature circle packing metric r av . 
27rx(X)^ > - Z(e,v)zLk(l)(* ~ He)) + 2n X (Fi) , VI : cf> ± I C V . 
Furthermore, if any of above four statements is true, the solution of combinatorial Calabi 



flow (2.1) converges exponentially fast to constant curvature circle packing metric v Q 



Proof: We have proved the long time existence of flow (2.1) in last section. Denote 



r(t), t £ [0, oo) as the solution of combinatorial Calabi flow (2.1) 



First we prove (1) => (3). If r(t) converges, i.e., r(+oo) = lim r(t) S R> exists, 
then both K(+oo) = lim K (t) E W and L(+oo) = lim L(t) exist too. Consider the 

t— >+oo t— >+oo 

combinatorial Calabi energy 



JV 



JV 



i=i 



i=l 



differentiate over t, we obtain 



N 



C'(t) = 2 ^ KiKi = 2K T K = -2K T L 2 K < 0, 

i=i 

which implies that C{t) is descending as time t increases. Moreover, C(t) is bounded by 
zero from below. Therefore C(+oo) = lim C(t) exists. C(+oo) exists too, since 

t->+oo 

C'(+oo) = lim (-2K T (t)L 2 (t)K(t)) = -2K T (+oo)L 2 (+oo)K(+oo). 

t— >+oo 
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The existence of C(+oo) and C'(+oo) implies C'(+oo) = 0, that is, 

K T (+oo)L 2 (+oo)K(+oo) = 0. 



Hence K (+oo) G Ker(L 2 ) = Ker(L). From Lemma 3.3 we know K av = K(-\-oo) G 
which implies that constant curvature metric r av exists. 

Next we prove (3) (1). Assuming there exists constant curvature circle packing metric 
r av , which implies K av G , we want to show r(t) converges, i.e. r(+oo) = lim r(t) G 

£—►+00 

H>0- We carry out the proof in three steps. 

Step 1: Denote Ai as the first eigenvalue of L, we claim that 

K T L 2 K> \\\\K - K av \\ 2 = X\C. (4.5) 

Select an orthogonal matrix P, such that P 7 LP = diag{0, Ai, • • • , Aat_i}. Write P = 
(eo, ei, • • • , ejv_i), where is the (i + l) th column of P. Then Leo = and Lei = A^, 
i = !,-■■ ,N —1, which implies e = -^(1, • • ■ , 1) T . Write K = P T K, then K\ = e^iT = 

vNk av . Hence we obtain 

K T L 2 K = \\K 2 + • • • + \ 2 N _ X K 2 N > \j(\\k\\ 2 - K 2 ) = \ 2 \\K - K av \\ 2 = \ 2 C. 

Thus we get the claim above. 

Step 2: we show that {r(t)\t G [0, oo)} CC E->q- Consider the combinatorial Ricci 
potential 

fU N 

/(")-/ ^2{Ki-k av )d Ui , uGR N , (4.6) 

J u av i=1 

where u av = La;p _1 (r at ,). This integral is well defined, because Y2i=i(^i ~ k av )dv,i is a 
closed differential form. /|^ is strictly convex (see the proof of Theorem B.2), hence 
V/l^ : 'Wa R N ,u •->• K — K av is injective, therefore u av is the unique critical point of 
and / is bounded from blow by zero (f(u) > f(u av ) = 0). Consider ip(t) = f(u(t)), 

then 

<p'(t) = (V/) T u(t) = (K — K av f(-LK) = —K T LK < 
hence ip(t) is descending as t increases and then < ip(t) < ty?(0) = f(u(0)), for any 
t G [0,+oo). Hence {u(t)\t G [0,+oo)} C [f\<%)~ ([0,(^(0)]). Because /|^ is proper 



by Theorem B.2, (f\ a ^) 1 ([0,(^(0)]) is a compact subset of % . Therefore {u(t) | t G 
[0, +oo)} CC or equivalently, 

{r(t)\te [0,oo)} cc ^ c R^ . 

Step 3: we show that r(t) converges exponentially fast to r av . Due to step 2, \f (t), the 
first eigenvalue of L 2 (i), has a uniform lower bound along the Calabi flow, i.e., Af (i) > 
A/2 > 0, where A is a positive constant. Using (4.5), we obtain 

C'(t) = -2K T L 2 K < -2\\(t)C < -AC. 
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Hence then we have 

C(t) < C{0)e- xt , (4.7) 

which implies that the curvature Ki(t) converges exponentially fast to k av . Moreover, we 
can deduce that r(t) and u(t) converge exponentially fast to r av and u av respectively by 
similar methods as in the proof of Lemma 4.1 in [12J. 

The equivalence between (2), (3) and (4) can be found in [8] (Thurston first claimed 
that (3) and (4) are equivalent in |26j). Thus we finished the proof. Q.E.D. 

Remark 3. We can prove {r(t) \ t € [0, T)} CC R>q similarly without assuming T = +oo. 
Using the classical extension theorem of solutions in ODE theory, we obtain T = +oo. 



This gives a new proof of Theorem 3.5 



Theorem 4.4. For (X, T, <£), consider the following user prescribed combinatorial Calabi 
flow 

u = L(K-K), (4.8) 
where K S H N is any prescribed curvature. For any initial circle packing metric r(0) £ 



R>0' the solution of (4.8) exists for t G [0, +oo). Moreover, the following three statements 



are mutually equivalent 



(1) 
(2) 
(3) 



The solution of (4.8) converges. 

The prescribed curvature K is admissible, i.e., K G . 

The solution of user prescribed Ricci flow u = K — K converges. 



Furthermore, if any of above three statements is true, the solution of flow (4.8) converges 
exponentially fast to f. 
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Appendix A. Dual structure and the proof of (3.2) 

A.l. Dual structure determined by circle patterns. Suppose Ci,Cj,Ck are closed 
disks centered at Vi, Vj and so that their radii are r^, rj and r^. They both intersect with 
each other at an angle supplementary to &(viVj), Q(vjVk), &(vkVi). Let £», Cj, Ck be the 
geodesic lines passing through the pairs of the intersection points of {Ck,Cj}, {Ck,Ci}, 
{Q, Cj}. These three lines £j, Cj, must intersect in a common point O (see Figure |3l). 
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EF={i,j} 



Figure 3. discrete dual-Laplacian 



Figure 4. dual structure 



Let i' be an intersection point of the circles Cj and Denote 



cos 9 



i'k 



ij | Aijk 



cos Of cos Q k - % 



sm< 



(A.l) 



where 9 % - k is the inner angle at Vj in the triangle {j, k,i'}. l*j\Aijk is the directed distance 
from O to H, the foot of perpendicular at edge ViVj. The directed distance from O to 
edge ViVj is positive (negative), when O is inside (outside) ZviVjVk- 

Think O as {i, j, k}* , the dual vertex to the triangle {i, j, k}. For two adjacent triangles 
k} and having a common edge {i, j}, {i, j}* (the dual edge to the edge {i, j}) 

is the line from {i,j, k}* to (see Figure]!]). Note that {i,j}* is perpendicular to 

{i, j}. Denote £?• as the length of {i, j}* , then I*- = l*j\Aijk + tfjlAijl- Using cosine law in 
triangle {i, j, k}, we have 

'' '''' i - . i [ i ^. i i . i i i ~~. i i . ' i « i — ^ r y . - ■ i 



dOf _ rj coB0 3 
dr^ 3 1 



A cos Of ■ cosOj l 



l* I 

''a | Aijk 



'J 



sin 



Qik 



I; 



Thus we obtain 



Bi 



7* I /* I /* 
,. _ Hj\Aijk HjlAijl _ ijj 

f)r- 3 r)r- 3 h I I 



Furthermore, 



N 



N 



A * = - E = - E ^/i = E ^(/i - /«) = E hf, - /<)> 



(A.2) 



(A.3) 



(A.4) 



'7 



i=i ' i=i i~* 

where / : V — > R is a function defined on vertices. 

The classical discrete Laplace operator "A" is often written as the following form ([101) 
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Figure 5. O is inside A ABC, such that H lies between B and C. 



where the weight ujij can be arbitrarily selected for different purpose. Thurston first 
studied t he du al structure of circle patterns in [26J . Bennett Chow and Feng Luo calculated 



formula (A. 2) in [SJ. Glickenstein wrote uiij = in [13]. The weight comes from dual 



structure of circle patterns. This is why we call (A.4) discrete dual-Laplacian. There are 



other types of discrete Laplacians, such as cotangent-Laplacian. See more in |15j , [16 
P2j, [ZU, and 122]- 



A.2. proof of (3.2). 



Lemma A.l. Consider a triangle AABC coming from a configuration of circle patterns 
with weight $ G [0, f]. Assuming O is inside AABC. OH is the altitude from O onto side 
BC. Assuming the point H lies between B and C (see Figure pi). Then |Oif| < \/3|BC|. 



\OH\ > V3\BC\ =>- < 



Proof: Suppose this lemma is false. Then \OH\ > ^/3\BC\, and 

\OH\ > V3\BH\ => IBOH < vr/6, 

< 

\OH\ > V3\CH\ =>- ICOH < vr/6. 



Hence 



and 



IBOC = IBOH + ICOH < tt/3, 



ZAOC = 2vr - + IBOC) > 2vr - (it + 

ZvlOS = 2vr - (ZylOC + ZSOC) > 2vr - (vr + 



2vr 
2tt 



In AAOB, we have > |vlO|, for Ai? faces to a much bigger angle. Thus we can select 
a unique point S between A and B so that \AS\ = \ AO\. using 



— = — < — 



IBAO vr 
"1 < 2' 
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we get 



lAOB-*,tAOB-**°» 
2' 2 



C 



7T 7T 

6' 2 



and 



/^c D tt-IBAO 

10SB = 7T G 



7T 






7T 


2vr" 


7T , 




c 






2' 


ff 2 


2' 


T 



Then we obtain 

which implies 
Similarly, we have 
Next we show 
and 



\OB\ 



\OB\ sin IB SO 1 

< , 



|AB|-L40| \BS\ Bin ISOB 



\AB\ - \AO\ > -\OB\. 
i ii 2 1 



2, 



- |AO| > -|OC|. 



r B > ^\OB\ 
rc>\\OC\. 



If fB < |OB| 3 we know < \OA\ by use of r\ — r\ = |OA| 2 — |05| 2 . We already know 
r A + r B > 1-^-^1 • Hence 

> |AB| - r A > \AB\ - \OA\ > -=\OB\. 

Thus we always have r# > i|0£?|, no matter r# > |Oi?| or < | Oi?| . Similarly, we have 
fc > ||OC|. Then it is easy to see 

r% + r 2 c > ^(\OB\ 2 + |OC| 2 ) > ~ • 2|0#| 2 = * |Otf| 2 > ^|£C| 2 > |5C| 2 , 
which contradicts with the fact 

\BC\- 
Thus the lemma is true. 



| + r 2 , + 2r B r c cos ^bc 1 >\lr%,+r 2 r . 



B'C 



Q.E.D. 



(Above Lemma belongs to Ruixiang Zhang and Chenjie Fan. 



Corollary A. 2. Given (A, T, <£), where A is a closed surface, T is a triangulation, $ 6 
[0,7r/2] is a weight. Then 



< 



drj ■' 



80. 



Proof: Thurston claimed that ^-rj = l*j\Aijk is positive ( [25] ) - Hence is inside 



the AviVjVk- The proof can be finished by Lemma A.l 
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Q.E.D. 



Appendix B. Combinatorial Ricci potential is proper 
Lemma B.l. Assuming function ip £ C 2 (E n ) satisfies 

(1) : tp is strictly convex to the downwards, and 

(2) : there exists at least one point p such that Vifj(p) = 0. 

Then lim ifj(x) = +oo. Moreover, tp is proper. 

x— >oa 

Proof: Set h(t) = inf ip(x), t > 0. Then h(t) is a strictly monotone increasing 

\x\=t 

function. We just need to prove h(t) —> +oo, as t — >■ +oo. The process is almost the same 
with Lemma B.l in [12], we omit the details. Notice that there is another method to prove 



Lemma B.l by remark 2.7 in [T]. Q.E.D. 



Theorem B.2. Given (X, T, <£), where X, T and $ are defined as before. Assuming there 

>JV 

L >o- 

N 



exists a constant curvature metric r av G E> . Consider the combinatorial Ricci potential 



/(«) =/E ( k * - k ™) dui 



where u av = Exp l (r av ). Then for arbitrary constant a, /L. is proper and 



lim f(u) = +oo. 



Proof: Write n = N — 1 for convenience, and denote a T = ^t^(1> • ■ "j 1) T as the unit 
normal vector perpendicular to the hyperplane Then we have f(u) = f(u + ta T ) for 
all t S E. Thus we only need to prove above theorem for a = 0. Write as % for short. 
The hyperplane % is determined by n linear independent variables. We claim that /|^ 
is strictly convex when considered as a function of n independent variables. 

Denote E™ = {( = (Cir - ' i(n) T £ R^ICjv = 0)}. Select an orthogonal linear trans- 
formation £/ : R N -)■ E^, such that g/a T = (0, • • • , 0, 1) T and, moreover, s/ transforms 
% to E? one by one. The matrix form of under standard orthogonal basis of E^ is 



denoted by A. Set c/ = o «e/ 1 

5(Cl,--- ,Cn)=/(V(Cl,--- ,Cn,0) 



T 



We want to show that Hess(g) is a n x n positive definite matrix. Partition the matrix A 
into two blocks 

" S ~ 



A 
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where S is a n x N matrix. By calculation we get V/ = K — K av , Hess(f) 



L, Vg = SVf = S(K — K av ) and Hess(g) = SLS T . Because Aa T 
have a = (0, • • • ,0, \)A = a^. Using La T = 0, we obtain 



d(K u ...,K N ) 
9(ui,...,njv) 



(0,--- ,0,1) 



wc 



ALA 



L [S 



SLS T SLa T ~ 




'SLS T 


" 


aLS T aLa T 











On one hand, SLS T is positive semi-definite, since ALA T is positive semi-definite. On the 
other hand, rank(SLS T ) = rank(ALA T ) = rank(L) = n. Therefore Hess(g) = SLS T 
is positive definite. Moreover, / (= go s^\^) is strictly convex when considered as a 
function of N — 1 variables. Thus we get the claim above. The proof can be finished by 
using Lemma B.l 

Q.E.D. 
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